Topologically Linked Polymers are Anyon Systems 
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We consider the statistical mechanics of a system of topologically linked polymers, such as for 
instance a dense solution of polymer rings. If the possible topological states of the system are 
distinguished using the Gauss linking number as a topological invariant, the partition function of 
an ensemble of TV closed polymers coincides with the 2N point function of a field theory containing 
a set of iV complex replica fields and Abelian Chern-Simons fields. Thanks to this mapping to 
field theories, some quantitative predictions on the behavior of topologically entangled polymers 
have been obtained by exploiting perturbative techniques. In order to go beyond perturbation 
theory, a connection between polymers and anyons is established here. It is shown in this way that 
the topological forces which maintain two polymers in a given topological configuration have both 
attractive and repulsive components. When these opposite components reach a sort of equilibrium, 
the system finds itself in a self-dual point similar to that which, in the Landau-Ginzburg model 
for superconductors, corresponds to the transition from type I to type II superconductivity. The 
significance of self-duality in polymer physics is illustrated considering the example of the so-called 
4 — plat configurations, which are of interest in the biochemistry of DNA processes like replication, 
transcription and recombination. The case of static vortex solutions of the Euler-Lagrange equations 
is discussed. 



In this letter we consider the statistical mechanics of 
a system of topologically linked polymers Q, such as for 
instance a dense solution of polymer rings (catenanes). 
If the possible topological states of the system are dis- 
tinguished using the Gauss linking number as a topolog- 
ical invariant, it has been shown in Q that the partition 
function of an ensemble of N closed polymers coincides 
with the 2N point function of a field theory containing 
a set of N complex replica fields and Abelian Chern- 
Simons (C-S) fields. Thanks to this mapping to field 
theories, some quantitative predictions on the behavior 
of topologically entangled polymers have been obtained 
by exploiting perturbative methods, see Ref. j3| and ref- 
erences therein. If one wishes to go beyond perturbation 
theory, however, one is faced with the problem of field 
theories which contain multi-component scalar fields and 
gauge fields interacting together. Clearly, it is not easy 
to study theories of this kind even in the simplest case 
N = 2. The idea which we wish to pursue here in or- 
der to circumvent such difficulties is to establish a link 
between polymers and anyons Q|. Anyon field theories 
have been in fact intensively investigated. It is known 
for instance that their classical equations of motion ad- 
mit vortex solutions Q . In a special region of the space of 
parameters, called the Bogoml'nyi self-dual point |fg, the 
attractive and repulsive forces between vortices vanish. 
This point corresponds in the Landau-Ginzburg model 
for superconductors to the transition from type I to type 
II superconductivity. 

The first obstacle in this program is that anyons are 
living in 2 + 1 dimensions, while polymer trajectories are 
intrinsecally defined in three dimensions. For this rea- 
son, one coordinate must be singled out and regarded 
as time. Accordingly, a point of the three-dimensional 
Euclidean space R 3 will be identified by a three vec- 



tor £ M = (r,i), where [i = 1,2,3, r = (x l ,x 2 ) and 
x-i = t plays the role of time. We need also to intro- 
duce the two-dimensional completely antisymmetric ten- 
sor e lJ , i — 1,2 and its three dimensional counterpart 
e A " yp , ix,v,p — 1,2,3. These tensors are uniquely de- 
fined by the following conventions: e 12 = 1 and e 123 = 1. 
We use middle latin letters as two-dimensional space in- 
dices and middle greek letters as three dimensional in- 
dices. Sum over repeated space indices will be every- 
where understood. The trajectories T a of the polymers, 
a = 1, . . . , N, will be treated as continuous curves. Usu- 
ally, these trajectories are parametrized by means of their 
arc-length a a : T a = {{r a (a a ),t a (a a ))\ < a a < L a }, L a 
being the total length of the trajectory. To make a con- 
nection with anyons, however, it is convenient to use the 
time t as a parameter. Of course, the coordinate t is 
not a good parameter if we do not introduce a suitable 
of sectioning for the curves r a . To this purpose, let us 
notice that, with respect to the t— direction, each curve 
T a has maxima and minima defined by the condition: 



dt a {a a ) _ 



0. The number of minima, s a , is equal to the 
number of maxima. Let us pick up a given point of mini- 
mum on F a and call it r a ,i. Starting from r a ,\ and going 
along the curve after choosing arbitrarily a direction, one 
will encounter successively a point of maximum, that will 
be called T a ,2, a point of minimum r Q) 3 and so on. In this 
way we obtain a set of 2s a points r a ,i, • ■ ■ , T a ^ Sa ■ Now it 
is possible to split T a into s a open paths r 0) i,... ,r ajSa 
which connect pairwise contiguous points of maxima and 
minima: 



Ta.I < t < T a j + i 

r a J a (T a j a+ i) — r a j a+ i(T a j a+ i) 
r a .i a (T a ,i a ) = r 0) / a _i(T /J 



(1) 
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where I a is a cyclic index such that I a € {1, ...,2s Q } 
and I Q + 2s a = I a . The points r a j a (T a j A ) are considered 
as fixed points. This sectioning procedure is explicitly 
illustrated in the case s = 3 by Fig.;Q] By construction, 



at the points of maxima and minima T a j a in order to 
reconstruct the closed curve r a . 

The polymer action A can be divided into two contri- 
butions: 




FIG. 1: Sectioning procedure for a 2s — plat F a with s = 3. 

there is a one-to-one correspondence between the points 
of the t— axis in the interval [r a ,i, t ,z+i] and the points 
of the trajectories r a j. As a consequence, it is possi- 
ble to parametrize the trajectories T a j a by means of the 
variable t. 

We are now ready to construct the partition function 
of a system of topologically entangled polymers. For the 
sake of clarity, it is convenient to use some of the nomen- 
cleature of knot theory. In particular, from the mathe- 
matical point of view, each closed trajectory T a can be 
regarded as a knot, while an ensemble of two or more 
trajectories defines a link. Moreover, with some slight 
abuse of language, a knot or a link with 2s points of 
maxima and minima will be called a 2s— plat. The con- 
cept of plats arises quite naturally in the processes of 
replication, recombination and transcription of DNA 

To start with, let us write the partition function for a 
link with the topology of a 25— plat, where S = J2a=i Sa - 
The integers s a are now kept fixed, so that the link is 
composed by knots T a which are 2s a — plats. In order 
to avoid a proliferation of indices, we will restrict our- 
selves to a system of two linked polymers, thus setting 
N = 2. Actually, this is not a big limitation, because 
the two-polymer problem contains all the ingredients of 
the N— polymer problem. Putting i — \Z^T, the desired 
polymer partition function is given by: 



Z(X) 



boundary 
conditions 



Vv aJa (t)e~ A (2) 



The boundary conditions on the trajectories T a ,i a are 
those of Eq. These boundary conditions are required 
by the fact that the open paths T a j a , with a fixed and 
I a = l,...,2s a , must join together in a suitable way 



A A free A 



The free part A free is: 
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free 



2 2s a 

EE 

o=l J„=l 



dt(-l 



top 



(3) 



9al a 



(*) 



(II 



(4) 

The factors (— l) Ia 1 appearing in Eq. iQJ are necessary 
to make Af ree positive definite. The parameters g a i a are 
related to the lengths L a Ia of the open chains r a j a as 
follows. Let us discretize T a Ia considering it as a random 
chains with n segments. If Eq. <@J would be the action of 
a set of free two-dimensional chains, in the limit of large 
values of n and keeping fixed the ratio n/g a j a , one would 
obtain the well known result 0]: L a j a ~ n/g a ,i a - Here 
things are a little bit more complicated, because one has 
to take into account the third dimension, which in turn 
is also the variable which parametrizes the curves. As a 
consequence, the above relation governing the length of 
the trajectories is replaced by: 



T 2 



\ T a,I a 



T aJa 



9 a, I a 



\Ta,I a + l — T a ,I a 



(5) 



The second term in Eq. {2J) is: 
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top 
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i=i 



dt 



dt 



B M (n,j(t),t) 



(6) 



In the above equation = (B,Bq) and C M = (G,Co) 
are two Abelian C-S vector fields. We have put B 3 = B n 
and C3 = Co in order to stress the role of time of the 
third spatial coordinate, k and A are real parameters. 
k coincides with the C-S coupling constant (see below). 
The spurious gauge degrees of freedom have been elimi- 
nated using the gauge condition (Coulomb gauge): 



V B = V C = 



(7) 



so that the gauge fixed C-S action Scs looks as follows: 

S cs = ^J d 3 i >V W'.. + Coe^diBj] (8) 

The C-S fields act as auxiliary fields which impose topo- 
logical constraints on the closed trajectories Ti and IV 
To show this, one has to integrate out these fields from 
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the partition function (2J and then to perform an inverse 
Fourier transformation of Z(A): 



(9) 



The details of this calculation have been already ex- 
plained in Ref. ;3j and will not be repeated here. The re- 
sult is the new partition function Z(m), which describes 
the fluctuations of two polymer rings constrained by the 
condition: 



in 



x(Ti,r 2 ) 



(10) 



where x^Fii^z) is the Gauss linking number and the 
topological number m must be an integer. A few com- 
ments are in order at this point. First of all, in 0] the 
topological constraint 11 Oil has been derived from the C- 
S path integral using the Lorentz gauge, while here the 
Coulomb gauge (JJJ has been chosen. Formally, things 
works exactly in the same way in both gauges. How- 
ever, one should keep in mind that the final expression 
of the Gauss linking number x(ri,r 2 ) obtained starting 
from the Lorentz gauge differs from the expression that 
one obtains starting from the Coulomb gauge. Of course, 
the two expressions are equivalent due to gauge invari- 
ance, but the way in which the winding of one trajectory 
around the other is counted, is realized in a very differ- 
ent way. Another important point is that our approach 
allows a more refined classification of links than the bare 
Gauss linking number. Indeed, two links which have the 
same Gauss linking number can be further distinguished 
by their plat-structure, i. e. by the number of maxima 
and minima of the knots which are composing the links. 
Finally, it is not difficult to add also the repulsive steric 
interactions between the monomers in the partition func- 
tion @. It is interesting to notice that, since we have 
chosen the coordinate x$ = t in order to parametrize the 
curves, the usual S— function potential which describes 
these forces takes the form: I-^r, t) — voS(r)S(t), vq being 
a temperature dependent constant. Due to the presence 
of the 5— function in time, the steric interactions become 
instantaneous in our case. In the following we will not 



discuss the steric interactions, in part because we prefer 
to concentrate on the pure topological interactions, in 
part because the inclusion of these repulsive forces com- 
plicates enormously the task of finding non-perturbative 
analytic solutions of the two-polymer problem under con- 
sideration. Thus, we will assume that in our system steric 
interactions are negligible. This approximation is valid 
for instance in polymer solutions which are very dilute 
or in which the temperature is near the so-called theta- 
point. The hypothesis of very low monomer densities is 
natural in the present context, because we are discussing 
the statistical properties of two isolated polymer rings, 
which are supposed to be very long and thin. 

To establish the connection with anyons, we need to 
pass in the partition function <|2jl from the statistical sum 
over the trajectories r a j a (t) to a path integral over fields. 
The details of the derivation of the partition function 
Z(X) in terms of fields in the general case and including 
also the steric interactions will be presented elsewhere. 
Here we will restrict ourselves to links which have the 
structure of 4— plats, putting henceforth s\ = s 2 = 1. 
This case is particularly interesting for biological appli- 
cations, since most of the links produced by in vitro en- 
zymology experiments are 4— plats [8j. Moreover, "small" 
links up to seven crossings are all 4 — plats. 

Introducing a set of n r complex replica fields: 



(ii) 



and a convenient notation for covariant derivatives: 
D(±A,B) = V±iAB D(±^,C) = V±z^-C (12) 

the 4— plat partition function can be expressed as follows: 

. 2 2 

Z 4 - P iat{\) = lim / V(fields) J[ ld J] 2 ,je- s ^^ 
n_> J i=i j=i 

(13) 

The action Sn-piat contains the free field action and the 
topological terms: 



Tl,2 



Si-plat 



T l,l 

T2,2 



dt dxi *J ^#1,1 



dt dx < ** 2 ,8^2,1 



2 l 

D(-A,B)*i,i +*i 2 St*i,2 



1 12 



D(A,B)*i, : 



D(- — ,C)* 2l i + tf£ 2 a t * 2 , 2 + D(— ,C)* 2 ,2 

2n ' Ag 2 2 27r 



(14) 



dx = dx\dx2 is the infinitesimal volume element in the two-dimensional space and dt = -Jr. Let us note that in 
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the partition function il'M we have already integrated out 
the time components of the Chern-Simons fields Bo, Co, 
which play the role of Lagrange multipliers and impose 
the following constraints: 



B = 2 



K 





2 






#2,1 




#2,2 







2 






#1,1 




#1,2 


1 



e(T 2 ,2-t)9(t-T 2 ,l)(l5) 
e(n t2 - t)9(t - Ti,i) (16) 



In the above equation B and C are the magnetic fields 
associated to the vector potentials B and C respectively: 



B = d x B 2 - d 2 B 1 = <"n,H, 
C = d 1 C 2 -d 2 C 1 =e l] d l C J 



(17) 
(18) 



and 9(t) is the Heaviside theta function 9(t) = if t < 
and 9{t) = 1 if t > 0. The symbol V(fields) in Eq. (O 
is a shorthand for the field integration measure: 

„ 2s a 

V(fields) = / II V$* aI V$ aJa (19) 

/a = l 

Finally, the operators O a ,i a are given by: 

O a J„ = "0a,/ a ( r a,/a( r a,/a + l), r aJa + l) 

^l(r a , /a (r a , /a ),T Qi/ J (20) 

Eqs. I| 1 3f) — (I20fl establish the desired connection between 
polymers and anyons. In fact, the action itTHl . together 
with the constraints lTL5l) - ljT6jl . coincides formally with 
the action of a multi-layered anyon system. The number 
of layers is 2(si + 82) = 4, while the fields z i^a,i a 



propagate anyonic particles with spin n r and the statis- 
tics of boson. One difference from anyonic field theories is 
due to the fact that, in the case of real particles, one con- 
siders the evolution of the whole system from an initial 
time T\ and a final time T2. Here, instead, time intervals 
measure the extensions of the trajectories r a ,/ a in the 
X3 direction. For this reason, in the partition function 
ltl3)l part of the system evolves within the time inter- 
val [71,1,1-1,2] and another part within the time interval 
[t2,i,T2,2]- This inhomogeneity in the time intervals and 
the use of the Coulomb gauge, which makes the interac- 
tions instantaneous, is the cause of the presence of the 
Heaviside 9— functions in Eqs. l|T5|l - irL6ll . In this sense, 
the complete analogy with anyons can be recovered only 
in the limit 



T2,l = Ti 



Tl,2 



T2,2 



T 2 



(21) 



Throughout the rest of this letter we will assume that 
the above condition is satisfied. This is to avoid technical 
complications in dealing with the 9— functions. 

The connection to the anyon problem suggest to apply 
to the action of Eq. 11411 the Bogomol'nyi transformations 
0. Using these transformation in a suitable way, it is 
possible to rewrite Si- V i a t in the following form: 



M-plat 



plat 



24, 



/ a t 



(22) 



Neglecting surface terms and introducing the notation 
D± = Di ± zL> 2 , where Di, i — 1,2, denotes the i— th 
space components of the covariant derivatives (1121 . one 
has that: 



Fi-plat 
1 

+ 



r-2 



dt I dx 



' 2 2 

o=l I a = l 



4S 



11 



£>+(-A,B)*i.i +- Ld+(A,B)#i, 2 

I 4512 I 



Mi 



21 



£>_(-— ,C)* a ,i 



4.9: 



22 



— ,C)* 2 .2 

Ztt 



(23) 



and 



Z 



A— plat 



dt I dx 



gn S12 



#2,1 



512 322 







2 






M 


#1,1 




#1,2 





(24) 



It is easy to recognize from Eq. l|23"jl that Ti-piat con- 
sists in the sum of the self-dual actions of four different 
families of anyons # Q ,j a ■ These families are coupled to- 
gether by the constraints ltl5|l - ltT6)l . The second contri- 
bution to Si.piati the interaction term J^-piat, contains 



instead Coulomb-like potentials, which attract or repel 
the monomers belonging to different trajectories T a j a . It 
is difficult to guess what is the influence of these Coulomb 
forces on the statistical behavior of the polymers. The 
problem is that their strength is proportional to A (see 
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Eq. 12411 1 and A is not a real coupling constant with a 
physical meaning. At the end, we are interested to elim- 
inate A and to express the polymer partition function in 
terms of the topological number m, as we did in Eq. 
m is in fact the true physical parameter which describes 
the topological state of the system. What is possible to 
conclude from Eqs. I|22ll - I|24|l is that topological forces 
have attractive and repulsive components, which inter- 
fere with the steric interactions. This result confirms at 
a non-perturbative level a previous result obtained using 
perturbative methods Q. Also in Ref. Q it has been 
shown in the limit in which fluctuations in the particle 
densities are relatively small that the C-S fields generate 
effective Coulomb interactions in anyon models. 

Remarkably, our approach reveals that the two- 
polymer system under consideration has a self-dual point. 
This point is reached when the parameters g a i a are equal: 



911 — 9l2 — 921 — 922 



(25) 



If the above relations are satisfied, it is easy to see that 
li-piat = 0, so that the only surviving terms in the ac- 
tion S^-piat are those coming from the self-dual contri- 
butions in T^-piat- We note that the existence of the 
self-dual point J25II is not submitted to the assumption 
of Eq. lj2lj. Indeed, let us suppose that the points of 
maxima and minima T a ,i a of the trajectories r a / a are 
not aligned as specified by Eq. 12H . The only difference 
in the expression of the interaction term 2i- p i a t is that 
now the values of T\ and T 2 , which give the boundaries 
of the integration over time, should be replaced by the 
following ones: 



Ti = max[ri,i,r 2 ,i] 



T 2 = min[ri i2 ,r 2 ,2] (26) 



The choice 112611 of integration limits in I^-piat is dictated 
by the presence of the 9— functions of Heaviside in the 
constraints ltl5ll - lfT6jl . Physically, Eq. ll26ll is motivated 
by the fact that, in the Coulomb gauge, all interactions 
become instantaneous. Therefore, if maxima and minima 
are not aligned, it is easy to see that the trajectories may 
only interact in the time interval [Ti, T 2 ], where T\ and T 2 
are given by Eq. ffify. Substituting the new integration 
limits (jSH in Eq. (EH, it is clear that, if Eq. (|25j) is 
satisfied, the effective Coulomb interactions in I^ p i at 
vanish identically, so that the action Si- P i a ts coincides 
once again with the pure self-dual part Ti- V i a t- Thus, 
to have self-duality we only need that the effects of the 
steric interactions are negligible, such as for instance in 
the case of very low monomer concentration or at the 
theta-point. 

The self-duality condition 12511 may be physically in- 
terpreted as the requirement that the trajectories T a j a 
must be homogeneous, in the sense that they should have 
equal persistence lengths. In this case, the attractive and 
repulsive forces described by the term ln-piat counterbal- 
ance themselves in a sort of equilibrium which establishes 



the self-dual regime. This does not mean however that 
attactions and repulsions between monomers disappear, 
because Coulomb-like forces still hide in the self-dual part 
of the action T^ v i at . 

We would like now to exploit further the analogy of the 
polymer problem with anyons field theories and search for 
possible non-trivial classical solutions which minimize the 
polymer free energy Si- P i a t- Of course, one should keep 
in mind that, even in the case of the simplest model of 
anyons, vortex solutions outside the self-dual point have 
been found only by means of numerical methods. Also 
at the self-dual point, there is no clear and detailed un- 
derstanding of the dynamics of C-S vortices [lQj. Here 
the situation is further complicated by the presence of 
replica fields and multiple families of anyons which are 
mixed together due to the constraints I|15J1 - H16J1 . In view 
of these limitations, it seems reasonable to restrict our- 
selves to the self-duality regime 1 2 o f) and to static solu- 
tions, which satisfy the conditions dttp^ Ia — for every 
value of a = 1, 2, I = 1, 2 and w = 1, . . . , n r . As a conse- 
quence, from now on we will consider the pure self-dual 
action: 



Si— pl a t 



£>+(-A,B)*i,i + £>+(A,B)*i, 2 



D-{-— ,C)# a ,i + D-(—, C)* 2 ,2 

Z7T Z7T 



(27) 



Of course, static solutions do not fit very well with the 
physical boundary conditions which should be imposed to 
the classical fields at the points of maxima and minima Ti 
and T 2 . For this reason, we are supposing here implicitly 
that polymers are very long and that the interval T 2 — T\ 
is large. In this way the relevant contributions to the free 
energy Si- P i a t come from instants t which are far enough 
from the extrema located at t = Ti and t = T 2 . At these 
intermediate points the trajectories fluctuate in such a 
way that, at each fixed instant t, they visit more or less 
with the same frequency the same locations of the two 
dimensional space spanned by the vectors r a ,i a (t). 

From the action I|27J1 one finds the following Euler- 
Lagrange equations of motion: 



D+(-A,B)Vr,i = D+(A,B)^, 2 = 



(28) 



D - ("5? C ) ^,1 = D - C ) M.2 = (29) 

To these equations one should also add the constraints 
(tl5ll - lfT6jl which determine the transverse components of 
the fields B,C and the Coulomb gauge fixing @. In 
total, after separating the real and imaginary parts in 
Eqs. il28fl - l|29[l . we obtain 12 x n r equations, which com- 
pletely specify the classical field configurations. It seems 
natural to try replica symmetric solutions of the kind: 



CI 



(30) 
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where a = 1,2, I a = 1,2 and w = 1, . . . ,n r . Moreover, 
we switch to polar coordinates in order to express the 
complex fields: 

tf«,/ a =e*-.'.pal (31) 
The consistency of Eqs. H28H - I|29[I requires that: 

0a, 1 = —<t>a,2 Pa,l = — — (32) 

Pa, 2 

Here ci and c 2 are arbitrary constants different from zero. 
Also the densities p a ,i a must not vanish. There are other 
ways to make Eqs. I|28|I - II2H|I consistent, but they lead ei- 
ther to unphysical or to trivial solutions. After imposing 
the consistency conditions (l32f) in Eqs. (I28ll - lj29jl . there 
remain still four independent relations, which can be used 
to fix the components of the vector fields B, C: 

\B % = fi^i H-ieyfl 1 ' In p M (33) 

K 1 

— Ci = di<j)2,i - -eij&hip^i (34) 

It is easy to prove that the Coulomb gauge requirement 
10 is satisfied only if the phases 4>\^i and 4> 2 , l are con- 
stant. Finally, inserting Eqs. I|33J1 - H34|1 in the constraints 
itlfill - lfTfijl . it is possible to determine the remaining de- 
grees of freedom pi ( i and p2,i- 

A In pi i = 4An r [ — p 2 ,i ) (35) 

\P2,1 J 



A\np 2:1 = 4An r (p 1A - (36) 

A further study of the above two equations requires nu- 
merical methods. 

Concluding, we have succeded in establishing a con- 
nection between topologically linked polymers and anyon 
field theories. We have limited ourselves to 4— plats, but 
it is possible to generalize our results to the case of any 
2s— plat. In our approach, links which belong to different 
2s— plat configurations are distinguished. This allows a 
better classification of the topological states of polymers 
than the bare Gauss linking number. In view of possi- 
ble applications to knot theory, it would be nice to ex- 
tend also to non-abelian C-S field theories our techniques 
based on the Coulomb gauge and on the sectioning pro- 
cedure of loops described above. 

The analogy between polymers and anyons has been 
exploited to investigate the statistical mechanics of 
4— plats, a particular class of links, which is relevant in 
biological applications. It turns out that, as anyon par- 
ticles in multi-layered systems, also 4— plats have a self- 
dual point. Within the hypothesis of Eq. ll2l)l . which 
demands that the points of minima T a ^i a +i are aligned 
at the same height t = Ti, while the points of maxima 



are aligned at the height t — T%, self-duality occurs when 
all open curves T a j a have the same length. Particular 
polymer configurations of this kind may exist in nature, 
for instance after the process of DNA replication. How- 
ever, we have seen that the requirements for the presence 
of a self-dual point are much more general: We only need 
that the effects of the steric interactions are negligible, 
such as for instance in the case of very low monomer 
concentration or at the theta-point. 

The existence of the self-dual point 12511 shows that the 
physics of interacting polymers is much richer than pre- 
viously expected. For instance, the classical equations of 
motion admit non-trivial solutions, which are the analogs 
of vortex-like solutions in anyon models. By analytical 
methods it is just possible to investigate static field con- 
figurations. These suffer the limitations mentioned after 
Eq. 1I27JI and cannot be defined if the requirements of 
Eqs. (I21|l are not fulfilled. In fact, if maxima and min- 
ima are localized at different heights on the t— axis, the 
action £>4_ p ; t becomes explicitly dependent on time due 
to the presence of the Heaviside 9— functions in the fields 
B and C. This time dependence is very mild. Its only 
effect is that the boundaries of integration over time in 
the various terms composing the action Si- P i a t do not 
coincide. Yet, this is sufficient to spoil the possibility of 
having purely static field configurations. For these rea- 
sons, any conclusion on the meaning of the static solu- 
tions in the case of polymers should be taken with some 
care. Having in mind these caveats, it is interesting to 
note that, in the solutions which we have studied here, 
it seems to prevail a repulsive force between couples of 
trajectories T a j a belonging to the same polymer. This 
is what Eq. <.°>2ll suggests, because the densities p a ,i are 
constrained to be inversely proportional to their counter- 
parts p a , 2 at each point x. 

It is worthing to point out that, despite their close 
resemblances, the mechanisms used to establish the self- 
dual regime in our polymer model and in usual anyon 
field theories are different. In anyon field theories, in fact, 
the existence of the self-dual point requires that the effec- 
tive Coulomb interactions, which arise after performing 
the Bogomol'nyi transformations, are canceled against 
true Coulomb interactions of charged particles |llj]. In 
the polymer model, instead, the effective Coulomb inter- 
actions, represented by the term Ii- P i a t, counterbalance 
themselves and the presence of other interactions is not 
needed. This self-balancing transition to self-duality may 
occur also in other physical problems where multi-layered 
systems of non-charged particles are relevant. Another 
application of our approach are polymer brushes. In that 
case, C-S field theories in the Coulomb gauge are able to 
take into account the effects related to the winding of 
neighboring polymer trajectories. 
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